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@ Interleaved digitizer array with calibrated sample timing. 



® A digitizer system includes M digitizers, each producing a separate waveform data sequence 
succession of instantaneous magnitudes of an input signal at sample times determined by a per od.c : ctock 
sign^smission of the dock signal to each digitizer is delayed by a corresponding ad,ustab e delay time so 
as to control the relative sample timing of the digitizers. To adjust sample t-m.ng. a sine wave s,gnal ,s apphed 
2 LTpu? signal to each digitizer such that me M digitizers produce M separate r^TZZTZZ 
Ssponse to id input signal and the M data sequences are interleaved and endowed o form> 
waveform data sequence. A first sequence of complex numbers representing a discrete Fbuner J*£*»"" *• 
single waveform data sequence is generated and then a second sequence of M complex numbers js itamrt 
torn elements corresponding to relative magnitude peaks of the first sequence. A third sequence of M I complex 
Ambers is generated^preLting an invense discrete Fourier transform of the second sequence and me phase 
angle of each number of me third sequence is computed and divided by me input s,gnal ^quencyto produce a 
sefof M numbers, each representing a timing error for a corresponding one of the M d.g.bzers. The time delay 
corresponding to each digitizer is men adjusted by me amount of the bming error. 
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INTERLEAVED DIGITIZER ARRAY WITH CALIBRATED SAMPLE TIMING 



Background of the Invention 



s The present invention relates in general to a high speed waveform digitizer system of the type which 
employs a parallel array of digitizers to sample and digitize a waveform at successive times, and in 
particular to a method and apparatus for measuring and adjusting the time interval between successive 
waveform samples. 

A typical waveform digitizer employs a sample and hold (S/H) circuit to sample and store a voltage 

70 representing the instantaneous magnitude of an input signal and an analog-to-digital (A/D) converter which 
produces digital output data proportional to the sample voltage, thereby representing the instantaneous 
magnitude of the input signal. The sampling rate of such a digitizer is usually limited by the rate at which 
the A/D converter can convert the sample voltage to a stable digital value. 

Digitizer systems comprising arrays of such digitizers have been developed which can sample and 

;5 digitize an input signal at a higher rate than any one digitizer can operate. The input signal is applied in 
parallel to each digitizer of the array. A single dock signal passes through a multiple tap delay line (or a 
succession of time delay circuits), each tap of the line providing the clock signal input to a successive one 
of the digitizers so that for each successive digitizer of the array, the clock signal is delayed by a longer 
time. Thus each digitizer samples and digitizes the input signal once during each cycle of the clock signal, 

20 but at a different relative time during the cycle. The sequences of output data produced by all of digitizers 
are collected and stored in an interleaved fashion to provide a single data sequence representative of the 
input signal. While the propagation delay of each A/D converter of each individual digitizer need only be 
slightly less than the period of the clock signal, the effective sampling rate of an array of M digitizers is M 
times the frequency of the clock si gnal. ^ . _ ■ 

^ Rir a clock signal of frequency VM, each section of the delay line should delay trie strobe signaTby 
1/f s so that samples are taken at evenly spaced intervals. The resulting data sequence is then equivalent to 
a sequence produced by a single digitizer operating at a sampling frequency of f s . The effective sampling 
rate of such a digitizer array (and therefore the sampling time resolution of the sampling system) can be 
increased by increasing the number M of digitizers in the array and the number of taps in the delay line, 

30 and decreasing the time delay between each tap of the delay lines. However, the accuracy with which tiie 
delay line can be adjusted to delay the clock signal places a limit on the sampling time resolution that can 
be obtained. As f s increases, the delay provided by each section of the delay line decreases, and small 
errors in the time delay provided by each section have an increasing effect on the relative accuracy wfth 
which each sample is timed. Differences in the response time of each digitizer, which can vary with ambient 

35 temperature or component aging, cause additional time delay errors. When the number of digitizers reaches 
a point at which the time delay errors become sub stantial in comparison to the nominal time delay of each 
delay line section, an increase in the number M of digitizers cannot substantially improve the accuracy with 
which the data sequence produced by the array characterizes the input signal. 

40 

Summary of the invention 



A digitizer system comprises an array of M similar digitizers, each producing a separate output data , 
45 sequence representing a succession of instantaneous magnitude samples of an input signal. A periodic v 
clock signal is applied to each digitizer to control sample timing, but the clock signal is delayed by an 
adjustable delay time prior to application to each digitizer. When the clock signal time delay associated with 
each delay circuit is properly adjusted, the separate output data sequences of the digitizers may be 
interleaved to form a single data sequence representing magnitudes of the input signal at times which are 
so spaced by intervals of l/f s seconds, where fj/M is the frequency of the clock signal. 

In accordance with the present invention, the adjustment is performed by applying a sine wave signal of 
known frequency f 0 as the input signal to the digitizer system and permitting the digitizer system to produce 
an interleaved waveform data sequenc comprising IM elements, where I is an integer, suitably 16. The IM 
elements of the data sequence are then windowed (weighted), suitably by a 4-term Blackman-Harris 
window, and the resulting windowed waveform data sequence is transformed by a discrete Fourier 

2 
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transform to produce a sequence of IM complex numbers representing a frequency spectrum of the 
windowed waveform sequence. A total of M elements spaced I elements apart and represents peaks of 
the frequency spectrum are then extracted from the frequency spectrum sequence to produce a first M- 
element sequence of complex numbers, which is then subjected to inverse Founer transform to produce a 
second M-element sequence of complex numbers. The phase angle of each number of the second 
sequence is determined and divided by the frequency of the input s.ne wave s.gnal. Each of the resulting M 
values indicates a timing error associated with a separate one of the dig.tizers. and the clock signal time 
delays are adjusted accordingly to eliminate the timing errors. 

The accuracy of the timing error measurement depends only on the accuracy with which the frequency 
of the input sine wave signal can be controlled, and not on the amplitude of the sine wave s.gnal. Since s.ne 
wave signal generators known in the art can produce sine wave signals of accurately controllable frequency, 
the timing errors can be accurately measured. 

It is accordingly an object of the invention to provide a method for accurately measuring sampling clocx 
sianal time delays in a digitizer system utilizing a digitizer array. 

It is another object of the invention to provide a high frequency array-type digitizer system in which 
sample timing is accurately controlled. , * • 

The subject matter of the present invention is particularly pointed out and distinctly claimed .n the 
concluding portion of this specification. However, both the organization and method of operation of the 
invention, together with further advantages and objects thereof, may best be understood by reference to the 
following description take in connection with accompanying drawings wherein like reference characters refer 
to like elements. 

Drawings 

FIG. 1 is a block diagram of a digitizer system utilizing a digitizer array; 

FIG. 2 is a time domain plot of a sine wave input signal to the digitizer array of FIG. 1 showing 
nominal^ample-timing;- 



30 FIG. 3 is a plot of the frequency spectrum of a sine wave; 

FIG. 4 is a plot of the frequency spectrum of an infinite data sequence produced by un.formly 

sampling and digitizing a sine wave; 

FIG. 5 is a plot of the frequency spectrum of an infinite data sequence produced by nonuniformly 

sampling and digitizing a sine wave; 
35 FIG 6 is a plot of the frequency spectrum of a rectangular function; 

FIGS. 7 and 8 are plots of the frequency spectrum of a data sequence produced by wmdow.ng a 
truncated, non-untformly sampled sinewave; and 

FIG. 9 is a flow chart for programming the computer of FIG. 1 in accordance with the present 

invention. 
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Detailed Description 

FIG 1 shows a block diagram of a digitizer system 10 adapted to dig.tize an input signal with a nominal 
sampling frequency f s . The digitizer system comprises an array of M conventional digitizers 12. each 
including a sample and hold (S/H) circuit 14 for caputring a voltage sample of an .nput signal connected to 
the digitizer through a switch 17 controlled by a computer 22. The S.H circuit 14 captures the voltage 
sample when strobed by a pulse from a strobe signal. An analog-to-digital (AD) converter 16 converts the 
voltage sample to representative digital data, and an acquisition memory 18 stores the A.D converter output 
data. Each digitizer 12 also includes a time base circuit 20 for producing the strobe s.gnal applied to S,H 
circuit 14 in response to each pulse of an externally generated clock signal. Time base circuit 20 also 
addresses and write enables acquisition memory 18. the memory address being incremented between each 
write operation. The time base circuit 20 begins strobing SH circuit 14 in response to the first clock signal 
pulse following receipt of a start command from computer 22 transmitted by way of a control and data bus 
24 interconnecting each digitizer 12 to computer 22. In accordance with instructions prov.ded by computer 
22 via bus 24, each time base circuit 20 may stop acquiring data after a predetermmed number of data 
values have been stored in acquisition memory 18 or after the input s.gnal reaches a predetermined trigger 
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level. Digitizers capable of performing the functions of digitizers 12 are often utilized in digital oscilloscopes 
for digitizing input signals. • 

Thus each digitizer 12 is adapted to store a sequence of data in its acquisition memory 18. eacn 
element of the data sequence representing the magnitude of an input signal at a different time as 
determined by each pulse of an input clock signal. The clock signal applied to each digitizer 12 is denved 
from a master clock signal of frequency (</M produced by a clock signal generator 15. The master clock 
signal passes through a series of unit delay stages 26, suitably comprising delay lines, each providing a 
nominal delay 1/I S . The master clock signal and the dock signal output of each delay stage are applied as 
inputs to separate adjustable delay circuits 28. and the output of each adjustable delay circuit 28 provides 
the clock input to a separate digitizer 12. The delay provided by each adjustable delay circuit 28 .s 
controlled by data supplied by computer 22. the amount of delay provided by each adjustable delay circuit 
being adjusted as described hereinbelow. 

FIG 2 illustrates how digitizing system 10 digitizes one cycle of an input signal, for example a sine 
wave, when the number M of digitizers 12 is 6 and the adjustable delay circuits are adjusted such that the 
delay between sampling by successive digitizers 12 is 1/f s . At time TO, a first pulse of the master clock has 
propagated through one adjustable delay circuit 28 to a first digitizer 12 (Digitizer 0) and has caused 
Digitizer 0 to sample the input signal. The data produced and stored in the acquisition memory of Digitizer 
0 in response to the voltage sample taken at time TO is represented in FIG. 2 by the length of arrow DO at 
time TO At time T1. the same master clock pulse has propagated through one delay circuit 26 and one 
adjustable delay circuit 28 to a second digitizer 12 (Digitizer 1) and has caused Digitizer 1 to sample the 
input signal. The resulting voltage sample is subsequently digitized and stored in the acquisition memory of 
Digitizer 1 as data value D1. As the first master clock pulse reaches each successive Digitizer 2-5. the 
digitizer samples the input signal and stores a digital representation 02-D5 of the sample. At time T6. the 
second pulse of the master clock reaches Digitizer 0 and it acquires and stores another data sample DO. At 
25 times T7-T9. the second master clock pulse reaches Digitizers 1-3 and they acquire and store other data 
samples D1-D3. 

Thus it may be seen that each pulse of the master clock causes each digitizer 12 to sample the input 
signal once, but at a different time. With the delay circuits 28 properly adjusted, the period between each 
sample taken by each digitizer is M/f 5 . but the period between samples acquired by successive digitizers of 

so the array is lrf,. At the end of a data acquisition cycle, computer 22 of FIG. 1 acquires the data stored in 
the acquisition memory of each digitizer 12 by way of bus 24 and sequentially stores the data in an internal 
memory the data being ordered in an interleaved fashion according to the relative time that each data value 
was acquired. Thus computer 22 assembles and stores a waveform data sequence representing the input 
signal which sequence is identical to a waveform data sequence that would be produced by a single 

35 conventional digitizer 12 if it could operate at a sampling frequency of f s . 

It is seen that an array of M digitizers, each operable at a maximum frequency of f^M. cooperate to 
digitize an input signal with a sampling frequency of f«. However, in order for digitizer system 10 to operate 
properly, the time interval between samples acquired by successive digitizers must be 1/f„ Adjustable 
delay circuits 28 are needed because digitizers 12 have response times which can vary with ambient 

40 temperature and com ponent aging. In addition, inaccuracies in delays provided by delay stages 26 and 
variation in lengths of various conductors or other components in the clock signal paths can also cause 
errors in sample timing. The adjustable delay circuits 28 permit the intervals between sampling by 
successive digitizer circuits to be properly adjusted to account for such sources of error. 

The adjustment is made by first determining the error in each time interval. To measure the error in 
accordance with the present invention, a sine wave signal generator 19 is connected to an alternate input 
terminal of switch 17 as shown in FIG. 1. and . switch 17 is operated so as to apply the output of signal 
generator 19 as an input signal to the digitizer system 10. Signal generator 19 produces a sine wave signal 
of a known, accurately controlled, frequency f„. Digitizer system 10 is then operated for I cycles of the 
master clock such that each of the M digitizers 12 produces and stores a waveform data sequence 

so comprising I data elements. The quantity I is an integer, suitably 16. Computer 22 then obtains and stores 
the M data sequences in an interleaved fashion to produce a single waveform data sequence compnsing IM 
elements. Computer 22 then windows (weights) the IM elements of the data sequence, suitably by a 4-term 
Blackman-Harris window, and the resulting IM element windowed sequence is transformed by a discrete 
Fourier transform to produce a sequence of IM complex numbers of magnitude representative of the 

55 frequency spectrum of the windowed sequence. Every Ith element of the frequency spectrum sequence is 
then extracted to produce a first M-element sequence of complex num bers. The first M-element sequence 
is then subjected to inverse Fourier transformation to produce a second sequence of M complex numbers. 
The phase angle of each number of the second sequence is determined and divided by the radian 
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frequency of the input sine wave signal. Each of the resulting M values indicates a re.at.ve om.ng error 
S wil a separate one of the digitizers, and computer 22 a,ters the control data supp.,ed to each 
Se dS drcuft a to change the delay provided by the delay circu,. by the amount of the assooated 
timing error $ by ^ ^ Fouf , er , form of , 

sineTave signal 1 of frequency f.. The spectrum includes spectral lines represented by verfca. arrows at 
LquencL - S the absolute value of the magnitude of a s,gna. component of frequences md.cated by 
Z o zonta. portion of a spectral line is represented by its length. The discrete Fourier 
q ( x) is an infinite sequence of complex numbers, the absolute amplitudes of a port.cn of which are plotted 
„ J RG^ The plot of FIG. 4 represents the convolution of the frequency spectrum of Fig. 3 w,th the 
freouencv spectrum of a periodic pulse signal of frequency f, However, the d.g.t.z.ng system 10 of FIG 1 
Suces 2 "dam sequence, rather man an infinite data sequence, and when the adjustable de.ay 
*£Z 28 are mp^ adjusted, the interval between successive samples ,s non-un,form. Consequen«y 
SXuen?y spectrum of a data sequence produced by digitizing system 10 by digrtumg an .nput s.ne 

' 5 ^TSJSE VZZgSL 10 digits an input sign, g<, it ^ uces a 

data sequence S = (g(to). git,). 0(h) Oft* 9<W.>. 1 comprising M interleaved subsequences So. *. 

.... and S M ., as follows: 



20 



25 



30 



s 0 = tg(t 0 ), g(t M ), g(t 2M ), • • *1 
s 1 = [glt^, g(t M+1 )» g(t 2 M+i ) ' * * #1 



S M _i s tg^M-l 1 ' 9 (t 2M-l>' 9( t 3M-l ) » #1 

35 it i S clear that the m-th subsequence S m is obtained by uniformly sampling the value g(t + U) at the rate 
?,MT Tr \ = 1/Mf! To construct the sequence S. one might insert (M-1, zeroes between samples ,n all 

f^^^^'^ o.O, ] and men shift the subsequence S . m posrtions to the 

- ¥F. m m2£«* o-i -i where 2 " is *• unit de,ay operator - 3,111 fma,ly sum 

up all the subsequences to obtain the original sequence, 

M-1 

S = £ V" m ' U) 
m=0 

The digital spectrum. G(w). of S can then be represented by summaAon of those of 5 m z"\ for m = 0 to M- 

50 1. 
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M-l GO 

G(w) = (1/MT) Z Z G a [w-k(2tr/MT)] 

m=0 k=-» 



j[v-k(2»/MT)]t B g-jrowT, {2) 



e 

10 Let r m be the ratio of mT-tm to the average sampling period T, i.e., let 
tn, = mT - r m T. (3) 
Then equation (2) may be rewritten as 

G(w) » (1/T) Z (1/M) Ve-^. w -M2./MT]r n T 
^s— od m=0 

e -jkn.(2ir/M) G a [w-k(2tr/MT) ] • (4) 

Equations (2) and (4) are two general representations of the digital spectrum of a non-uniformly sampled 
signal. 

For a sinusoidal input signal g(t) = & w o l of frequency f 0 where w 0 = 2*f 0 , the Fourier transform is 
given by 

G a (w) = 2y6(w -w 0 ) (5) _ _.. 



30: 



where 5 is the delta function. 

Substituting equation (5) into equation (4). we have 



M-l « 

G(w) = (1/MT) Z Z 2ir6[w-W 0 -k(2ir/MT) J 

35 TO-0 k«-» 

e -jr m 2i,f 0 /f s e -jkm(2ir/M) (6) 

40 

where f s is the average sampling frequency 1/T. Let us define a sequence A(k), k = 0. 1. 2 M-1. M. * • * as 

follows. 

A(k) -V [l/M e-3 r m 2 *V f s] e -W 2 «'«>' < 7) 
m=0 

^ Then we can rewrite equation (6) as 

G(w) (1/T) I A(k)2ir 5[w-w Q -k(2»/MT)]. (8) 
k=— oo 

55 

Equations (7) and (8) are the complete spectraJ representation of a non-uniformly sampled sinisoidal signal. 
FIG. 5 shows a portion of the frequency spectrum of an infinite data sequence acquired by an M = 6 
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digitizer array in response to a sine wave input signal of frequency f. = f,4M. From equaflo <7> .t.s seen 
She sequence A(k) is periodic on k with the period M. hence the spectrum G(w, given by equatior (8V s 
oeriod'c on w with the period equal to 2*T = 2*f s , i.e.. the average sampl.ng radian frequency. 
SeLre one period of the spectrum comprises M line spectra uniformly spaced on the frequency axis 
s rSTSiSiS spectra, lines separated by the amount of f,M. The mam signal componen ,s ocated 
at 7 a.d I with magnitude proportional to |A(0)|, while the m-th harmonic component ,s located at ^ ^ 
and Ss maqrtSe proportional to |A(m)|. It is also noted that the relative strength among afferent 
and has magnitude ' P^ 0 "' , A , k .,, s Li e xaaaer ated in RG. 5 to indicate that all A(k)'s are. .n general. 

* 22. r" 3SS Uon m. m ■ «•■*■ < 81 " 

» shc^d be noted tt>at .0, each value 0. k tee-a is only one valoe .1 w to, «o,ch tee del* lonofon , non- 
zero in other words. G(w) is equal to f s A(k) for one particular value of k. 

ft is not« Tfram equation (7) that the finite sequence [A(k). k-0. 1 M-1] ,s the d.screte Fourier 

transform of the sequence 

is [(1/M)e-' r m 2,rf o'' , s],m=0. 1 M-1]. 

Given an IM-element sequence of complex numbers representing the spectrum G(w) of^uation (8). one 
could extract elements thereof to form an M-element sequence [w^k), k-0. 1....M-11- For example, tor 
^STg the spectrum shown in FIG. 5. one might extract sequence elements representing 
qSw^Z .... G(5V6 ♦ f 0 ) to produce sequence w s t A(0). A(1). A(2). A(3). A(4). A(5)l. By taking the 
discrete inverse Fourier transform of this sequence, we obtain the sequence 

wje-i^irfo/f, . m=0. 1 M-1]. 

If we then compute the phase angle of each element of this sequence we have the sequence: 

[2ftJAm=0, 1. .... M-1]. 
-•Oivldlng-each-^toment ofthe abov-sequ««eby2-f*«H»-no«n9Hhat-T. «~i.f < ..we*btain.the. M queoce: - 

[.nJ.n^O. 1 M-1]. 

Since r m T represents the time de.ay error of the Mth adjustab.e delay circuit 28 of RG. 1, this sequence 
indicates the amount by which each adjustable delay circuit must be adjusted. FIR , hv 

Thus one might determine the adjustment error of each adjustab.e de.ay arcu.t 28 o HG_ 1 I by 
anoMna a sine wave input signal of frequency f 0 to the digitizer system, taking the d.screte Founer 
KZ. of ZrSn 'sequel, extracting the sequence proportion, to A(k, ,o, r K -0 i 
me phase angle of each element of the sequence and dividing each phase angle by 2-V 
aTove described method assumes that the data sequence representing the input ^« °f '"finrte length 
whereas the data sequence produce by digitizer system 10 of FIG. 1 is of finite length. One may think of 
2 w^efom f data sequence produced by system 10 as representing sampled values of a waveform 
represent nTtheToduc* of a rectangu.ar .unction and the input sign,, trie rectangu.ar funct.cn hav.ng a un,t 
• a^Ttor a time interv, T and zero value before and after the interval. The frequency spectrum of the 
^tangujar Son. shown in FIG. 6. includes a major .obe of width 2<T centered about zero frequency and 

m Te Squ'ety 5 fpT^ St-" -a sequence produced by digitizer system ,0 would 
therefore not appear as a series of discrete spectr, lines as shown in FIG. 5. but wouW appear as the 
»^on ofle spectrum of F.G. 5 with the spectrum of F.G. 6. By windowing j^htrng terms < ) a 
reTanqTa? function suitably with a four-term Blackman-Harris window, the ampl.tude of the side lobes of 
me sMCfrum of FIG. 6 can be substantia.* reduced. The B.ackman-Harris window and other w.ndow.ng 
EcEfU described in the artide "On the Use o. Windows for Harmon.c Ana.ysis wrth ^.*=r*e 
TraSorm". by Frederic J. Harris, in Proceedings of the IEEE. Vo.. 66. No. 1 . January 1974. wh.ch arfc.e is 

^Ts^raCu^spectrum produced by taking the discrete Fourier transform of a waveform 
data Squence acquired by digitizer system 10 of FIG. 1 when M=6. and where the input sign, is a s.ne 
wave oTfZency \ - V* and the resu.ting sequence has been windowed by a B.ackman-Harns w.ndow^ 
TnTabsolute value of the amp.itude function A(f) has been p.otted on a .oganthm.c scale along the > vertc, 
L*Ts and frequency f is plotted along the horizont, axis. Major .obe pairs (shown ,n dotted hne) are centered 
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about the frequency of each spectrum line of FIG. 5 and have an amplitude equivalent to the amplitude of 
the corresponding spectrum line of FIG. 5. (Minor lobes are so small as to be negligible and am therefore 
not shown in FIG. 7.) The 4-term Blackman-Harris window is chosen for its low side-lobe and easy 
implementation. However, this window has a somewhat larger main lobe than other popular windows. In 
s order to avoid the interference from neighboring main lobes the input test frequency must be at least (1/4M) 
f s away from (m/M)f Sl where m = 1, 2. ... M. in one direction and at least (a'4M)f, away from those 
frequencies in the other direction. Hence at least 16M data points are needed. The absolute magnitude of 
each data sequence element of the transformed, windowed sequence is shown in FIG. 8 as a heavy line 
comprising the sum of the functions shown in dotted line. 
io Since the waveform data sequence produced by system 10 of FIG. 1 includes 16M elements 
representing the input waveform g(t). a 16M element sequence of complex numbers representing a. 
frequency spectrum of the input signal is produced by windowing the waveform data sequence with the 
Blackman-Harris window and taking the discrete Fourier transform of the result With M = 6. the sequence 
includes 16x6 = 93 complex numbers, each number representing a separate value G(f) where f=0, f*64. 

J5 2f y 64 63V64. The fourth number of the sequence of 96 complex numbers represents 6(4**96) = G(fo) 

and is equal to w^O). Similarly the (16k + 4)th number of the sequence represents G(4kfs<96) which is 
equal to w s A(k). Thus an M-element sequence wJA(m), m=0. 1. .... M-1] may be extracted from the 96- 
etement frequency spectrum sequence by choosing the 4th element of the frequency spectrum sequence 
and every 16th element thereafter. The timing error sequence [r m T m=0, 1. .... M-1] can th?n be obtained 
20 by determining the phase angle of each element w s A(m) as described hereinabove. 

If the waveform data sequence were not windowed, the minor lobes associated with each major lobe 
shown in FIG. 7 would not be negligible and would adversely affect the values of A(k) obtained and 
therefore the computed values of r m T would be less accurate. The windowing function could have been 
chosen to provide greater than or less than 50% overlap of the major lobes as shown in FIG. 7. However 
with greater than 50% overlap, the data values of G(f) which correspond to peak values of the lobes and 
which represent values of A(k). would be more strongly affected by adjacent major lobes. With less than 
50% overlap spectral resolution is decreased. 
- ^e frequency4 0 -of-the test signal- should-not^ 

avoid aliasing. In addition, in order to ensure that spectral peaks are evenly spaced and appear at 
frequencies corresponding to data points of the transformed, windowed data sequence G(f). f 0 should be 
selected to satisfy the following relation: 

fo = ( Vn) - V4M (n =M . 2. 3. ...) [10] 

In the example of FIG. 8. n = 12 which gives f c = V24. However, since the timing errors r m T are computed 
from the phase angles (27rr m f A), to maximize the resolution with which the phase angles are determined, it 
is desirable to make f 0 as large as possible so that the phase angles will be as large as possible. To do so 
it is necessary to choose n as small as possible. If n = 1 . then from equation [1 0], f G = 23f ff 24. But to avoid 
aliasing f Q must be not greater than tj2. Choosing n = 2. from equation [10] gives f G = 11 V24. which is 
less than half f s . Therefore the choice of n = 2 in equation 10 gives the sine wave signal frequency providing 
the best resolution in the determination of timing errors r m T. 

FIG. 8 shows the frequency spectrum of the windowed Fourier transform of a data sequence acquired 
from an M = 6 digitizer array acting on an input sinewave of frequency f Q ■ 11 ^24. The frequency 
spectrum is similar to that shown in FIG. 7. but has been shifted to the left by 10f y 24. As in the case of the 
45 frequency spectrum of FIG. 7. the sequence w^m) may be obtained by extracting the 4th frequency 
spectrum sequence number and every 16th number thereafter, but the resulting sequence is wJAW, A{5). 
A(0). A(1), A(2). A(3)] which must be reordered to obtain the sequence w s A(m), m = 0 to 5. 

FIG. 9 is a flow chart for programming computer 22 of FIG. 1 to measure and adjust the time delay 
circuits 28 of digitizer 10. Starting in step 40. delay circuits 28 of FIG. 1 are adjusted to provide an average 
so delay midway between minimum and maximum. Switch 17 is set to apply the sine wave output of sinewave 
signal' generator 19 as the input signal to digitizer system 10 (step 42). Each digitizer 12 is set to digitize 16 
samples of the input signal, and data acquisition is initiated (step 44). After the digitizers 12 have each 
acquired and stored sequences of 16 data values in their acquisition memories, computer 22 reads the data 
sequences and interleaves them in memory to form a single 96 element waveform data sequence (step 46) 
55 which is then windowed by a 4-term Blackman-Harris window (step 48). The discrete Fourier transform of 
the windowed sequence is then computed (step 50) to form a first data sequence, and a second data 
sequence of M elements is formed by extracting the 4th element of the first sequence and every 16th 
element thereafter (step 52). A third data sequence comprising M complex numbers is then computed (step 
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54) by taking the inverse discrete Fourier transform of the second data sequence. Each element of the third 
data sequence may then be normalized to the element corresponding to k = 0. i.e.. the element having a 
phase angle from which the error associated with Digitizer 0 is computed (step 56). The phase angle of 

• each element of the third data sequence is computed (step 58) and then divided by the radian frequency of 

5 the sine wave input signal (step 60) to produce the set of M time delay error values. Computer 22 then 
adjusts each time delay circuit 28 in accordance with the error values. It should be noted that the 
normalizing step (step 56), may be eliminated. However, the normalizing step obviates the need to aajust 
the delay circuit 28 associated with Digitizer 0. 

While a preferred embodiment of the present invention has been shown and. descnbed. it will be 

io apparent to those skilled in the art that many changes and modifications may be made without departing 
from the invention in its broader aspects. The appended claims are therefore intended to cover all such 

- changes and modifications as fall within the true spirit and scope of the invention. 



iff Claims 

1 For a digitizer system comprising an array of M digitizers each responsive to an input signal and a 
periodic clock signal of known frequency, where M is an integer greater than 1. and means for delaying 
transmission of the clock signal to each digitizer by a separate delay time corresponding to said each 

20 digitizer, each digitizer comprising means for producing a separate waveform data sequence in response to 
the clock signal, each data element of the separate waveform data sequence representing an instantaneous 
magnitude of the input signal during each period of said clock signal, a method for determining a timing 
error in the delay time corresponding to each digitizer, the method comprising the steps of: 

applying a sine wave signal of known frequency as the input signal to each digitizer such that the M 

25 digitizers produce M separate waveform data sequences in response to said input signal and said clock 
signal; 

generating a single waveform data sequence in accordance with a combination of the M separate 
waveform data sequences, the single waveform data sequence being representative of said input signal; 
generating a- first-sequence-of -complex numbers representing-a Jrequency .spectrum _of.the..single„ 

30 waveform data sequence; 

generating a second sequence of M complex numbers by extracting M elements representing relative 

magnitude peaks of the first sequence; 

generating a third sequence of M complex numbers representing an inverse discrete Founer transform 

of the second sequence; 

generating a set of M phase angle numbers, each phase angle number representing a phase angle 
associated with a separate complex number of the third sequence and corresponding to a separate one of 
said digitizers; and 

determining a sample timing error for each digitizer in accordance with the corresponding phase angle. 

2. The method in accordance with claim 1. wherein the step of generating a single waveform data 
sequence in accordance with a combination of the M separate waveform data sequences comprises the 
substeps of: 

interleaving elements of the M data sequences to form an interleaved data sequence: and 
windowing with a windowing function the interleaved data sequence to form said single waveform data 
sequence. 

3. The method in accordance with claim 1 wherein said sine wave signal has a frequency not greater 
than one half the frequency of said clock signal. 

4. The method in accordance with claim 3 wherein said clock signal has a frequency f, and said sine 
wave signal has a frequency f 0 satisfying the expression 
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so f 0 = (Vn) - (V4M) 

where n is selected from among the set of all integers greater than 1 . 

5. The method in accordance with claim 3 wherein said clock signal has a frequency U and said sine 
wave signal has a frequency f 0 satisfying the expression 
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6 The method in accordance with claim 1 wherein the step of determining a sample timing error for 
each digitizer in accordance with the corresponding phase angle comprises the substep of dividing the 
corresponding phase angle by an amount proportional to the frequency of said input signal. 

7. For a digitizer system comprising an array of M digitizers each responsive to an input signal and a 

5 periodic clock signal of known frequency f s , where M is an integer greater than 1. and means for delaying 
transmission of the clock signal to each digitizer by a separate delay time corresponding to said each 
digitizer, each digitizer comprising means for producing a separate waveform data sequence in response to 
the clock signal, each data element of the separate waveform data sequence representing an instantaneous 
magnitude of the input signal during each period of said clock signal, a method for determining a timing 

to error in the delay time corresponding to each digitizer, the method comprising the steps of: 

applying a sine wave signal of known frequency f Q as the input signal to each digitizer such, that the M 
digitizers produce M separate waveform data sequences in response to said input signal and said clock 
signal, where f 0 satisfies the expression f 0 = (f* n) - (f ff '4M) where n is selected from among the set of all 
integers greater than 1 ; ■ 

75 interleaving elements of the M data sequences to form a single interleaved data sequence 

representative of said input signal; 

windowing the interleaved data sequence to form a windowed waveform data sequence; 

generating a first sequence of complex numbers representing a discrete Fourier transform of the single 
waveform data sequence; 

20 generating a second sequence of M complex numbers comprising M elements of said first sequence 

corresponding to relative magnitude peaks of the first sequence; 

generating a third sequence of M complex numbers representing an inverse discrete Fourier transform 
of the second sequence; 

generating a set of M phase angle numbers, each phase angle number representing a phase angle 
25 associated with a separate complex number of the third sequence and corresponding to a separate one of 
said digitizers; and 

determining a sample timing error for each digitizer in accordance with the corresponding phase angle 
- number. ........ 

a The method in accordance with claim 7 wherein said sine wave signal has a frequency f G of 
30 frequency satisfying the expression 

f 0 = (V2) - (V4M). 

9. The method in accordance with claim 7 wherein the step of determining a sample timing error ror 
each digitizer in accordance with the corresponding phase angle number comprises the substep of dividing 

35 the corresponding phase angle number by an amount proportional to the frequency of said input signal. 

10. A self-calibrating digitizer system comprising 

an array of M digitizers each responsive to an input signal and a periodic clock signal of known 
frequency, where M is an integer greater than 1. each digitizer comprising means for producing a separate 
waveform data sequence in response to the clock signal, each data element of the separate waveform data 
40 sequence representing an instantaneous magnitude of the input signal during each period of said clock 
signal; 

means for delaying transmission of the clock signal to each digitizer by an adjustable delay time 
corresponding to each digitizer; 

means for generating a sine wave signal of known frequency for application as the input signal to each 

45 digitizer such that the M digitizers produce M separate waveform data sequences in response to said input 
signal and said clock signal; and 

means responsive to said M separate waveform data sequences for generating a single waveform data 
sequence in accordance with a combination of said M separate waveform data sequences, the single 
waveform data sequence being representative of said input signal, for generating a first sequence of 

so complex numbers representing a frequency spectrum of the single waveform data sequence, for generating 
a second sequence of complex numbers by extracting M elements representing relative magnitude peaks 
of the first sequence, for generating a third sequence of M complex numbers in accordance with an inverse 
discrete Fourier transform of the second sequence, for generating a set of M phase angle numbers, each 
phase angle number representing a phase angle associated with a separate complex number of the third 

55 sequence and corresponding to a separate one of said digitizers, for determining a timing error correspond- 
ing to each digitizer in accordance with the corresponding phase angle number, and for adjusting the 
adjustable delay time corresponding to each digitizer in accordance with the timing error corresponding to 
each digitizer. 
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11. A self-calibrating digitizer system comprising: 

an array of M digitizers each responsive to an input signal and a penod.c clock signal of known 
frequency f. where M is an integer greater than 1. each digitizer compnsing means tor producing a 
serrate waveform data sequence in response to the clock signal, each data element of the separate 
waveform data sequence representing an instananeous magnitude of the input signal during each penod of 

said clock signal; . . „ ..„. 

means for delaying transmission of the clock signal to each digitizer by a corresponding delay time 
means for generating a sine wave signal of known frequency f 0 as the input signal to each digitizer 
such that the M digitizers produce M separate waveform data sequences in response to said input signal 
and said clock signal, where f 0 satisfies the expression f 8 = (f r n) - (f r 4M) with n is selected from among 
the set of all inteqers greater than 2; and 

means for interleaving elements of the M data sequences to form a single interleaved data seouence 
representative of said input signal, for windowing the interleaves data sequence to form a wnoowea 
waveform data sequence, for generating a first sequence of complex numbers representing a discrete 
Fourier transform of the single waveform data sequence, for generating a second sequence ol M complex 
numbers by extracting M elements representing relative magnitude peaks of the first sequence, for 
generating a third sequence of M complex numbers representing an inverse discrete Founer transform of 
L second sequence, for generating a set of M phase angle numbers, each phase angle . number 
representing a phase angle associated with a separate comptex number of the h,rd sequence .and 
conesponding to a separate one of said digitizers, for determining a timing error for each digitizer in 
accordance with the corresponding phase angle number, and for adjusting the adiustable delay time 
conesponding to each digitizer in accordance with the corresponding timing error. 

1 1 The serf-calibrating digitizer system in accordance with claim 1 1 wherein sad sine wave signal has 
a frequency f 0 satisfying the expression 

fo = (tVn) - (V4M). 
where n is an integer greater than 1 . 
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0 Interleaved digitizer array with calibrated sample timing. 



© A digitizer system includes M digitizers, each 
producing a separate waveform data sequence re- 
presenting a succession of instantaneous magni- 
tudes of an input signal at sample times determined 
by a periodic clock signal. Transmission of the clock 
signal to each digitizer is delayed by a correspond- 
ing adjustable delay time so as to control the relative 
sample timing of the digitizers. To adjust sample 
timing, a sine wave signal is applied as the input 
signal to each digitizer such that the M digitizers 
produce M separate waveform data sequences in 
response to said input signal and the M data se- 
quences are interleaved and windowed to form a 
single! waveform data sequence. A first sequence of 
complex numbers representing a discrete Fourier 
CO transform of the single waveform data sequence is 
^generated and then a second sequence of M com- 
qgptex numbers is formed from elements correspond- 
r"ing to relative magnitude peaks of the first sequence. 

third sequence of M complex numbers is gen- 
COerated representing an inverse discrete Fourier 
0> transform of the second sequence and the phase 
^ angle of each number of the third sequence is 
O computed and divided by the input signal frequency 
^to produce a set of M numbers, each representing a 
qj timing error for a corresponding one of the M 
digitizers. The time delay corresponding to each 
digitizer is then adjusted by the amount of the timing 
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